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Abstract 

We revisit the question whether the running-mass inflation model allows the formation 
of Primordial Black Holes (PBHs) that are sufficiently long-lived to serve as candidates 
for Dark Matter. We incorporate recent cosmological data, including the WMAP 7-year 
results. Moreover, we include "the running of the running" of the spectral index of the 
power spectrum, as well as the renormalization group "running of the running" of the 
inflaton mass term. Our analysis indicates that formation of sufficiently heavy, and hence 
long-lived, PBHs still remains possible in this scenario. As a by-product, we show that the 
additional term in the inflaton potential still does not allow significant negative running 
of the spectral index. 
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1 Introduction 



The Cosmic Microwave Background (CMB) is very smooth. Full-sky observations allow to 
expand the measured CMB temperature in spherical harmonics Y% m . One can then determine 
the size of the anisotropies as a function of £, with larger i corresponding to smaller angles, 
and hence smaller length scales. Current CMB observation probe this power spectrum down 
to (comoving) length scales of about one Mpc. These observations imply very small primordial 
density perturbations at such large length scales, characterized by the power V-r, c — 10~ 9 . 

However, it is possible that the primordial density perturbations become much larger at 
smaller length scales, beyond the range probed by cosmological observations. Indeed, it is con- 
ceivable that these perturbations are so large that overdense regions collapse to form Primordial 
Black Holes (PBHs) just after the end of inflation [UEJE]. They are called "primordial" since 
they do not originate from the gravitational collapse of burnt-out stars; they could thus have 
any mass, including masses well below or well above stellar masses. Here we assume that the 
fluctuations which sourced the PBHs are also generated during inflation, specifically towards the 
end of inflation, well after the length scales probed by conventional cosmological observations 
exited the horizon. 

There are various constraints on PBH formation. For example, the density of roughly stellar 
mass black holes has to satisfy limits from searches for microlensing. Very light black holes could 
have evaporated (via Hawking radiation) in the epoch of Big Bang nucleosynthesis, altering the 
predicted isotope abundances. These and other constraints have recently been compiled in [4]. 
They can be translated into upper limits on the amplitude of the power spectrum at the length 
scales relevant for PBH formation, typically Vn c < 10~ 2 — 1CT 1 with some scale dependence. 

Clearly the power spectrum has to change dramatically towards the end of inflation for PBH 
formation to occur. In the framework of standard slow-roll inflation, this implies that the slow- 
roll parameters, and hence the inflaton potential, also should show large variations. One simple, 
and yet well motivated, model that can feature large variations of the slow-roll parameters is 
the running- mass model [5J E], a type of inflationary model which emerges naturally in the 
context of supersymmetric extensions of the Standard Model. The model is of the single-field 
type, but nevertheless it has relatively strong scale-dependence (running) of the spectral index. 
Previous studies of this model [7] showed that it can indeed accommodate sufficiently large 
positive running of the spectral index to allow for PBH formation. However, the recent 7-year 
WMAP data [S] prefer negative running of the spectral index. 

In this paper we therefore revisit the running mass model. We expand on refs.[7] by in- 
cluding the running of the running of the spectral index, which is only weakly constrained 
by existing CMB data. Similarly, we include the running of the running of the inflaton mass 
parameter in the potential. We show that for some (small) range of parameters, this model 
can still accommodate sufficiently large density perturbations near the end of inflation to allow 
for the formation of PBHs with mass larger than 10 15 g, which are sufficiently long-lived to 
be candidates for Cold Dark Matter (CDM). As a by-product we show that even with the 
additional term in the potential, the model does not allow for significant negative running of 
the spectral index. 

The remainder of this paper is organized as follows: In section 2 we briefly review the for- 
malism of PBH formation. In section 3 we discuss the running-mass inflation model, including 
the running of the running of the mass. We compute the slow-roll parameters which allow us 
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to determine the spectral index, its running, and the running of the running. In section 4 we 
perform a numerical analysis where we compute the spectral index at PBH scales exactly, by 
numerically solving the equation of motion of the inflaton field. Finally, we conclude in section 



2 Formation of Primordial Black Holes 

PBHs may have formed during the very early universe, and if so can have observational 
implications at the present epoch, e.g. from effects of their Hawking evaporation [1] for masses 
<• 10 15 g, or by contributing to the present "cold" dark matter density if they are more massive 
than 10 15 g [3]. (These PBHs would certainly be massive enough to be dynamically "cold"0). 

The traditional treatment of PBH formation is based on the Press-Schechter formalism 
[TO] used widely in large-scale structure studies. Here the density field is smoothed on a scale 
R(M). In the case at hand, R(M) is given by the mass enclosed inside radius R when R crossed 
the horizon. The probability of PBH formation is then estimated by simply integrating the 
probability distribution P(S; R) over the range of perturbations 5 which allow PBH formation: 
5 t h < 5 < 5 cut , where the upper limit arises since very large perturbations would correspond to 
separate closed 'baby' universes [SHU]- We will show that in practice P(5; R) is such a rapidly 
decreasing function of 5 above S t h that the upper cutoff is not important. The threshold 
density is taken as 5 t h > w i where w = p/p is the equation of state parameter describing the 
epoch during which PBH formation is supposed to have occurred [3]. Here we take w = 1/3, 
characteristic for the radiation dominated epoch which should have started soon after the end 
of inflaton. However the correct value of the threshold 5th is quite uncertain. Niemeyer and 
Jedamzik [12] carried out numerical simulations of the collapse of the isolated regions and found 
the threshold for PBH formation to be 0.7. We will show that PBHs abundance is sensitive to 
the value of <5 t h- 

The fraction of the energy density of the Universe in which the density fluctuation exceeds 
the threshold for PBH formation when smoothed on scale R(M), 8{M) > S t h, which will hence 



Here P(S; M(R)) is the probability distribution function (PDF) of the linear density field 5 
smoothed on a scale R, and 7 is the fraction of the total energy within a sphere of radius R 
that ends up inside the PBH. 

1 It has been argued [9] that BH evaporation might leave a stable remnant with mass of order of the Planck 
mass, which could form CDM; we will not pursue this possibility here. 

2 Throughout we assume for simplicity that the PBH mass is a fixed fraction 7 of the horizon mass corre- 
sponding to the smoothing scale. This is not strictly true. In general the mass of PBHs is expected to depend 
on the amplitude, size and shape of the perturbations ]V2 \ [T3 ] . 

3 We follow ref. [14] in including a factor of two on the right-hand side; this factor is not very important for 
delineating the inflationary model parameters allowing significant PBH formation. Moreover, we set 5 cu t to 
infinity. 



5. 





(2.1) 
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For Gaussian fluctuations, the probability distribution of the smoothed density field is given 

mR) =v^r p (-^m)- (2 - 2) 



This PDF is thus uniquely determined by the variance as{R) of 5, which is given by 

a 2 5 (R) = J W 2 (kR)V s (k)^-. (2.3) 

In order to compute the variance, we therefore have to know the power spectrum of 5, Vs{k) = 
fc 3 /(27r 2 ) (|<5fc| 2 ), as well as the volume-normalized Fourier transform of the window function 
used to smooth 5, W{kR). 

It is not obvious what the correct smoothing function W(kR) is; a top-hat function has 
often been used in the past, but we prefer to use a Gaussian window function^, 

W{kR) = exp . (2.4) 



Finally, on comoving hypersurfaces there is a simple relation between the density perturba- 
tion S and the curvature perturbation 1Z C [16J: 



2(1 + w) ( k x 2 



The density and curvature perturbation power spectra are therefore related by 



4(1 + wf ( k 



•i 



p ' (M > = (Tin#UJ r «' (ky (2 ' 6) 

The mass fraction of the Universe that will collapse into PBHs can now be computed by 
inserting eqs.f l2.6p and f !2.4[) into eq. fl2.3p to determine the variance as function of R. This has 
to be used in eq.( l2.2p . which finally has to be inserted into eq.f l2.ip . Since we assume a Gaussian 
P{$) in eq.f l2.2p . the integral in eq. fl2.ip simply gives an error function. 

In order to complete this calculation one needs to relate the mass M to the comoving 
smoothing scale R. The number density of PBHs formed during the reheating phase just after 
the end of inflation will be greatly diluted by the reheating itself. We therefore only consider 
PBHs which form during the radiation dominated era after reheating is (more or less) complete. 
The initial PBHs mass Mpbh is related to the particle horizon mass M bjHI 

4:71 

M PBR = 7 M = — 1P H- 3 , (2.7) 

when the scale enters the horizon, R = (aH)~ l . Here the coefficient 7, which already appeared 
in eq. fl2.ll) . depends on the details of gravitational collapse. A simple analytical calculation 



4 This PDF is often written as P(6(R)). However, we think it is more transparent to consider P to be the 
PDF of S, which is just an integration variable in eq. (|2.1[) . Eq. (|2.2p shows that the functional form of P(S) 
depends on the parameter i?, which in turn depends on the horizon mass M. 

5 Bringmann et al. |15j argued that a top-hat window function predicts a larger PBH abundance. 

6 Throughout the paper we put c = h = ke = 1- 
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suggest that 7 ~ w 3 ^ 2 ~ 0.2 during the radiation era [3]. During radiation domination aH oc 
a -1 , and expansion at constant entropy gives p oc g* ^ 3 a~ 4 [TTJ (where g* is the number of 
relativistic degrees of freedom, and we have approximated the temperature and entropy degrees 
of freedom as equal). This implies that 

\ i/3 



where the subscript "eq" refers to quantities evaluated at matter-radiation equality. In the 
early Universe, the effective relativistic degree of freedom g* is expected to be of order 100, 
while g^eq = 3.36 and k eq = 0.07fi m /i 2 Mpc _1 (fi m /i 2 = 0.1334 [8]). The horizon mass at 
matter-radiation equality is given by 

M 7T-3 _ 4?r Prad,0 , s 

q ~~ ~^P^m H eq - -y fc3 I > l^- y J 



eq°eq 



where a e q = (1 + ^ eq ) = 3146 and (assuming three species of massless neutrinos) f2 r ad,o^ 
4.17 x 10 -5 . Then it is straightforward to show that 

R _ _. _i /M PRH \ 1/2 / o* \V6 



1 Mpc 



Note that M PB h oc R 2 , not oc i? 3 as one might naively have expected. Recall that M PB h is 
related to the horizon mass at the time when the comoving scale R again crossed into the 
horizon. Larger scales re-enter later, when the energy density was lower; this weakens the 
dependence of Mpbh on R. Moreover, the lightest black holes to form are those corresponding 
to a comoving scale that re-enters the horizon immediately after inflation^ 

The Gaussian window function in eq. fl2.3p strongly suppresses contributions with k > 1/R. 
At the same time, the factor k 4 in eq. f)2.6p suppresses contributions to the integral in eq. fl2.3p 
from small k. As a result, this integral is dominated by a limited range of k— values near 1/R. 
Over this limited range one can to good approximation assume a power-law primordial power 
spectrum with fixed power Vn c (k) = Vn c (k R )(k/ k^ 3 ^' 1 , with k R = 1/R. With this 
ansatz, the variance of the primordial density field at horizon crossing is given by 



2 _ 2(l + w) 2 , 



s{R) = (5 + 3wY Vnc{kR)T[{ns{R) + 3)/2] ' (2 ' n) 
for n s (R) > -3. 

The power Vn c is known accurately at CMB scales; for example, P^ c (/co) = (2.43 ±0.11) x 
10 -9 at the "COBE scale" ko = 0.002 Mpc^ 1 [8J. In order to relate this to the scales relevant 
for PBH formation, we parameterize the power spectrum as 

Vn c (k R ) = Vn^kn/kor^- 1 . (2.12) 



7 In fact, PBH formation might also occur on scales that never leave the horizon [T5]. We do not consider 
this contribution here. 

8 "Fixed" here means that ns does not depend on k; however, ns does depend on R, since a large range of 
values of R has to be considered for PBH formation of different masses. 
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It is important to distinguish between ns{R) and n(R) at this point. ris(R) describes the slope 
of the power spectrum at scales k ~ ftR — 1/R, whereas n(R) fixes the normalization of the 
spectrum at k^ 3> k . The two powers are identical if the spectral index is strictly constant, 
i.e. if neither n$ nor n depend on R. However, in this case CMB data imply [8] that n = n$ is 
close to unity. Eqs. fl2.lip and f 1 2 . X 2 j) then give a very small variance, leading to essentially no 
PBH formation. 

Significant PBH formation can therefore only occur in scenarios with running spectral index. 
We parameterize the scale dependence of n as [19]: 

n(R) = n s (k )-^a s \n(k R) + ^(3 s ln 2 (k R) + ... ; (2.13) 

recall that we are interested in R <C l/k , i.e. \ia(k R) < 0. The parameters as and j3s denote 
the running of the effective spectral index ns and the running of the running, respectively: 



n s (k ) 
a s (k ) 
0s(k o ) 



dliak 

dns 



k=ko 



dink 

d 2 ns 



d\n 2 k 



(2.14) 

k=ko 



Eq.f l2.13j) illustrates the difference between n(R) and ns(R). The latter has an expansion 
similar to eq.f l2.13p . but with the usual Taylor-expansion coefficients, 1 in front of as and 1/2 
in front of (3$- One therefore has 

n s (R) = n(R)-^a s \n(k R) + ^ s ^ 2 (k R) + ... . (2.15) 

Setting n s (k ) = 1 for simplicity, eq. f)2.15p implies n s (R) = 2n(R) — 1 for 0s = 0, and 
ns(R) = 3n(R) — 2 for as = 0. We will compute the variance o~(R), and hence the PBH 
fraction /, for these two relations; they represent extreme cases if neither as nor (3s is negative. 

The result of this calculation is shown in figure [TJ Here we have fixed 7 = 0.2, and show 
results for two choices of the threshold 5 t h and three choices of n(R). We see that scenarios 
where n(R) = 1.3 (or smaller) are safe in the SM, because there is no model-independent limit 
on / for Mpbh < 10 10 g [1]. As noted earlier, PBHs contributing to Dark Matter today must 
have Mpbh^IO 15 g; at this mass, they saturate the DM relic density if / ~ 5 x 10 19 Figured] 
shows that this requires n(R) ~ 1.37 (1.41) for <5 t h = 0.3 (0.7). The dependence on ns is much 
milder. 

Figure [T] also illustrates a serious problem that all scenarios that aim to explain the re- 
quired CDM density in terms of post-inflationary PBH formation face. We just saw that this 
can happen only if the spectral index n increases significantly between the scales probed by 
the CMB and other cosmological observations and the scale R ~ 10~ 9 pc relevant for the for- 
mation of 10 15 g PBHs. However, n must then decrease rapidly when going to slightly smaller 



9 Note that / describes the fraction of the energy density in PBHs at the time of their formation. Since 
they behave like matter at all times, their fractional contribution to the energy density increases during the 
radiation-dominated epoch, and stays essentially constant during the subsequent matter-dominated epoch. 
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Figure 1: Fraction of the energy density of the universe collapsing into PBHs as a function 
of the PBH mass, for three different values of n(R) and two different choices of the threshold 
5 t h = 0.3 (0.7) for the solid (dashed) curves. On the upper [lower] of two curves with equal 
pattern and color we have assumed ns(R) = 2n(R) — 1 [ns(R) = 3n(R) — 2]. 

length scales, since otherwise one would overproduce lighter PBHs. For example, successful Big 
Bang Nucleosynthesis requires [I] /(10 13 g) < 2 x 10 -20 , about 12 orders of magnitude below 
that predicted by keeping n(R) fixed at the value required for having 10 15 g PBHs as CDM 
candidates. We will come back to this point at the end of our paper. 

Another problem is that the expansion of n(R) in eq. fl2.13l) will generally only be accurate 
if | ln(fco-R)| is not too large. This is the case for cosmological observations, which probe scales 
> 1 Mpc. The expansion becomes questionable for the scales probed by PBH formation. 
For example, fixing k = 0.002 Mpc" 1 , eq. fl2TT0|) gives |ln(fc -R)| = 41.1 for M PBH = 10 15 g. 
Fortunately within the framework of a given inflationary scenario this second problem can be 
solved by computing n(R) and ns(R) exactly, rather than using the expansion (I2.13p . 

Observational bounds on % and as can e.g. be found in ref.JS]. By requiring that these 
parameters remain within the 2a ranges for all scales k between 10 -4 Mpc -1 and 10 Mpc -1 , 
we find 0s ~ 0.03. Here we used the error bars derived from the analysis of WM AP7 data plus 
data on baryonic acoustic oscillations (BAO), ignoring possible tensor modes (as appropriate 
for small-field models), and including independent measurements of the Hubble constant as 
a prior (i.e. we used the "WMAP7+BAO+ifo" data set). Here we use the pivot scale k = 
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0.0155 Mpc -1 where ns(k ) and as(k ) are uncorrelated. On the other hand, eg. (12. 131) shows 
that for as = we only need (3s(k ) ~ 0.0016 in order to generate sufficiently large density 
perturbations to allow formation of 10 15 g PBHs. Even if we set «s(fco) equal to its central value, 
a s{ko) = —0.022, we only need (3s(k ) ~ 0.0032. Including the running of the running of the 
spectral index thus easily allows to accommodate PBH formation in scenarios that reproduce 
all current cosmological observations at large scales. 

Of course, this kind of model-independent analysis does not show whether simple, reason- 
ably well-motivated inflationary models exist that can generate a sufficiently large fls- One 
scenario that quite naturally accommodates strong scale dependence of the spectral index is 
the running-mass inflation model [5j 16]. We now apply our results to this model. 



3 Running- Mass Inflation Model 

This model, proposed by Stewart [5], [6], exploits the observation that in field theory the 
parameters of the Lagrangian are scale dependent. This is true in particular for the mass of the 
scalar inflaton, which can thus be considered to be a "running" parameter!*"! The running of 
the mass parameter can be exploited to solve the "77— problem" of inflation in supergravity [20J. 
This problem arises because the vacuum energy driving inflation also breaks supersymmetry. 
In "generic" supergravity models the vacuum energy therefore gives a large (gravity-mediated) 
contribution to the inflaton mass, yielding |?7| ~ 1. However, this argument applies to the scale 
where SUSY breaking is felt, which should be close to the (reduced) Planck scale Mp = 2.4 x 10 18 
GeV. In the running mass model, renormalization group (RG) running of the inflaton mass 
reduces the inflaton mass, and hence \rj\, at scales where inflation actually happens. There are 
four types of model, depending on the sign of the squared inflaton mass at the Planck scale, 
and on whether or not that sign change between Mp and the scales characteristic for inflation 

The simplest running-mass model is based on the inflationary potential 

V^Vo + ^mK^ 2 , (3.1) 

where is a real scalar; in supersymmetry it could be the real or imaginary part of the scalar 
component of a chiral superfield. The natural size of \m 2 i (Mp)\ in supergravity is of order 
Vo/Mp. Even for this large value of m 2 ^ which gives \r)\ ~ 0(1), the potential will be dominated 
by the constant term for (ft 2 <C Mp; as mentioned above, running is supposed to reduce |m|| 
even more at lower <p 2 . 

The potential (13.11) would lead to eternal inflation. One possibility to end inflation is to 
implement the idea of hybrid inflation To that end, one introduces a real scalar "waterfall 
field" ip, and adds to the potential the termj"1 



V^^-J^ + fr 1 . (3.2) 



10 Note that the physical, or pole, mass of the inflaton is not "running"; however, at the quantum level the 
physical mass differs from the parameter appearing in eq. (|3.f I) even if <fi = 0. 

n The parameters of this potential will in general also be scale dependent; however, this is immaterial for our 
argument. 



s 



Here A and k are real couplings, and the coefficient of the ip term has been chosen such that 
Vi n { = Vcj, + has a minimum with (V) = if = 0, ip ^ 0. As long as 2 > \/^§r, i> 
remains frozen at the origin. Once 2 falls below this critical value, ip quickly approaches its 
final vacuum expectation value, given by (ip) 2 = \J"^r-i while quickly goes to zero, thereby 
"shutting off" inflation. However, in this paper we focus on the inflationary period itself; the 
evolution of perturbations at length scales that left the horizon a few e-folds before the end of 
inflation should not be affected by the details of how inflation is brought to an end@ 

During inflation, the potential is thus simply given by eq. (l3.ip . Here m 2 (0) is obtained by 
integrating an RG equation of the form 

dm 2 , 

SE* = A " (3 - 3) 

where f3 m is the (3— function of the inflaton mass parameter. If m 2 is a pure SUSY-breaking 
term, to one loop (3 m can be schematically written as [SJ [21] 

2C D 

m = am 2 + — — | Ay | 2 m\ , (3.4) 

7T I07T Z 

where the first term arises from the gauge interaction with coupling a and the second term 
from the Yukawa interaction Ay. C and D are positive numbers of order one, which depend 
on the representations of the fields coupling to 0, m is a gaugino mass parameter, while m 2 is 
the scalar SUSY breaking mass-squared of the scalar particles interacting with the inflaton via 
Yukawa interaction Ay. 

For successful inflation the running of m 2 , must be sufficiently strong to generate a local 
extremum of the potential for some nonvanishing field value, which we call 0*. The inflaton 
potential will obviously be flat near 0*, so that inflation usually occurs at field values not very 
far from 0*. We therefore expand m 2 (0) around = 0*. The potential we work with thus 
reads: 

V = V + l -m 2 A^ 2 + ~ c 2 In (± )+ ~ 9 <P In 2 ( 4~) . (3.5) 



Here c 



dm 2 



<iln< 



d 2 ml 

is given by the p— function of eq. (!3.4l) . and g = 



is given by 



d(ln0) 2 

the scale dependence of the parameters appearing in eq. (13.41) . In contrast to earlier analyses 

of this model [2TI |24"| 125] . we include the In 2 (0/0*) term in the potential. This is a two-loop 

correction, but it can be computed by "iterating" the one-loop correction^ Since the coefficient 

g of this term is of fourth order in couplings, one will naturally expect \g\ |c|. However, 

this need not be true if |c| "happens" to be suppressed by a cancellation in eq. (l3.4p . Including 

the second correction to wi 2 (0) seems natural given that we also expanded the running of the 

spectral index to second (quadratic) order. 

Recall that we had defined 0* to be a local extremum of V^, i.e. ^'(0*) = 0. This implies 
1 

[2T] m 2 = m?(0*) = — — c; this relation is not affected by the two-loop correction oc g. 



12 It has recently been pointed out that the waterfall phase might contribute to PBH formation [25] . 
13 There are also "genuine" two-loop corrections, which can not be obtained from a one-loop calculation, but 
they only affect the term linear in ln(</>/0*). They are thus formally included in our coefficient c. 
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In order to calculate the spectral parameters as and 0s defined in eqs.f l2.14l) . we need 
the first four slow-roll parameters, defined as [IE] 





Ml fV'\ 


e 


2 \V) 


V 


- M *v ' 


e 


- Mp V 2 ' 


a 3 


R V' 2 V"" 

= M i v 



(3.6) 

where primes denote derivatives with respect to 0. All these parameters are in general scale- 
dependent, i.e. they have to be evaluated at the value of 4> that the inflaton field had when 
the scale k crossed out of the horizon. The spectral parameters are related to these slow-roll 
parameters by [16J: 



n s = 1 - 6e + 2n , 
a s = -24e 2 + 16erj - 2£ 2 , 
0s = -192e 3 + 192e 2 r/ - 32erf 



24er + 2r]C + 2a 3 . 



(3.7) 



Slow-roll inflation requires |e|, \r]\ 1. Combining eqs.f l3.5p and ( 13. 6ft we see that we need 
Vo 3> c<f) 2 L, g<p 2 L 2 , where we have introduced the short-hand notation 



In 



(3.8) 



In other words, the inflaton potential has to be dominated by the constant term, as noted earlier. 
Eqs. (13.61) then imply two strong inequalities between (combinations of) slow-roll parameters: 



|e| 

1 



< \v\ 

« \e 



(3.9) 



The first relation means that ns — 1 is essentially determined by n. Similarly, both relations 
together imply that as is basically fixed by £ 2 , while only the last two terms in the expression 
for /3s are relevant; these two terms are generically of similar order of magnitude. Finally, the 
factors of V appearing in the denominators of eqs.f l3.6j) can be replaced by Vq. The spectral 
parameters are thus given by: 



n s - 1 



a s 



cMl 



Vo 



L + 1 + (L 2 + 3L + 1) 



-2 



cMl 
V 



L 



cMl 
~Vo~ 



l + ^(2L + 3) 



— (L 

2c V 



fc {L + V > 



9 



^(3L + 2) + ^(3L 2 + 5L + - 



(3.10) 



14 



The powers on ^ 2 and a 3 are purely by convention; in particular, £ 2 could be negative. 
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where L has been defined in eq. fl3.8p . Clearly the spectral index is not scale-invariant unless 
c and g are very close to zero. Note that Vq appears in eqs. fl3.10p only in the dimensionless 
combination cMp/Vo, while only appears via L, i.e. only the ratio 0/0* appears in these 
equations. 

In contrast, the absolute normalization of the power spectrum is given by (in slow-roll 
approximation) 

1 V 3 

Vn < = u^mIv^ • (3 ' n) 

This normalization is usually quoted at the "COBE scale" k$ = 0.002 Mpc -1 , where V-r, = 
2.43 x 10~ 9 |8J. Applying our potential (I3.5P to eq.f l3.lip . replacing V by Vo in the numerator, 
we see that Vk c not only depends on cM|/Vo and L, but also on the ratio Vb/(Mp0 2 ). We can 
thus always find parameters that give the correct normalization of the power spectrum, for all 
possible combinations of the spectral parameters. 

We want to find out whether the potential (13.51) can accommodate sufficient running of 
ns to allow PBH formation. There are strong observational constraints on «j and as- It is 
therefore preferable to use these physical quantities directly as inputs, rather than the model 
parameters cM|/Vq, L and g/c. To this end we rewrite the first eq.f l3.10p as: 



cM| n s - 1 



Vo 2(L + 1) + 9 - (L 2 + 3L + 1) 
Inserting this into the second eq. fl3.10p gives: 



(3.12) 



a s = _^l L [l + U2L + ,)] [l + *(L + l)] 

2 [ L + l + JL (L 2 + 3L + 1 y2 



We thus see that the running of the spectral index is "generically" of order (ns — l) 2 ; similarly, 
the running of the running can easily be seen to be cx (ns — l) 3 . This is true in nearly all 
inflationary scenarios that have a scale-dependent spectral index. 

Eq. fl3.13l) can be solved for g/c. Bringing the denominator to the left-hand side leads to a 
quadratic equation, which has two solutions. They can be written as: 
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(L + 1) (L 2 + 3L + 1) + r s L (1.5L + 2) ± LJr s [{L + l) 2 + l] + r 2 (0.5L + 1 



2c " (L 2 + 3L + l) 2 + r s L (2L + 3) (L + 1) 

where we have introduced the quantity 



(3.14) 



(n s ~ 1) 



2 



r s = ^ '- . (3.15) 

2a s 

Since g and c are real quantities, eq. (I3.14p only makes sense if the argument of the square root 
is non-negative. Note that the coefficients multiplying rs and r| inside the square root are 
both non-negative. This means that the model can in principle accommodate any non-negative 
value of rs- However, small negative values of rs cannot be realized. It is easy to see that the 
constraint on rs is weakest for L = 0. The argument of the square root is then positive if 



11 



2r s + r| > 0, which implies either r$ > or r$ < —2. Recalling the definition (13. 15ft we are 
thus led to the conclusion 



as > 



(3.16) 



this bound should hold on all scales, as long as the potential is described by eq. (I3.5D . Note that 
it is identical to the bound found in ref. (25], i.e. it is not affected by adding the term oc L 2 
to the inflaton potential. This is somewhat disappointing, since recent data indicate that as 
is negative at CMB scales. Even the generalized version of the running mass model therefore 
cannot reproduce the current 1 a range of ag. 

However, at the 2 a level significantly positive as values are still allowed. Let us therefore 
continue with our analysis, and search for combinations of parameters within the current 2 a 
range that might lead to significant PBH formation. Using eqs. (13.121) and (13.141) we can use 
ns(ko) — 1, as(k ) and L = ln(0 o /0*) as input parameters in the last eq. (13. 101) to evaluate 
Psiko). This can then be inserted into eq. (12. 131) to see how large the density perturbations at 
potential PBH scales are. 

This numerical analysis is most easily performed at the "pivot scale", where the errors on 
n s and as are essentially uncorrected; at values above (below) this scale, n s and as are 
correlated (anticorrelated) [26J. The pivot scale for the "WMAP7+BAO+iJo" data set we are 



using is ko = A: pivot = 0.0155 Mpc 1 jE]0 At this scale, observati 



ions give 



n s {k ) = 0.964 ±0.013; 
a s (k ) = -0.022 ±0.020. 



(3.17) 



We saw above that requiring the correct normalization of the power spectrum at CMB scales 
does not impose any constraint on the spectral parameters. However, the model also has to 
satisfy several consistency conditions. To begin with, it should provide a sufficient amount of 
inflation. In the slow-roll approximation, the number of e-folds of inflation following from the 
potential (13.51) is given by: 



AN(L) 



c Jl L> [I 
2 

0^1) 



In 



L 



1 + 



L 



In ■ 2c 



o 



(3.18) 



where we have introduced the dimensionless quantity 



2cMg 
Vn 



(3.19) 



recall that it can be traded for ns(k ) — 1 using eq. (l3.12j) . Moreover, L = ln(0/ '0*) can be 
related to the scale k through 

k(L) = k e ANiL) . (3.20) 



15 This is the smaller of two fc p i vo t values given in ref. [3- The difference between these two values is not 
important for our numerical analysis. 



12 



This can be inverted to give 



L(k)-L E Wd + £) (321) 

m - L °l + llL (l-E(k)) + l]> (3 ' 21) 

where we have introduced 

E(k) = (-) . (3.22) 

The problem is that the denominator in eq. (13.211) vanishes for some finite value of k. This 
defines an extremal value of AN: 



AiV cx = -, r In 

5(1 + £) 





2c V 








~J). 



(3.23) 



A negative value of AN ex is generally not problematic, since only a few e-folds of inflation 
have to have occurred before our pivot scale k crossed out of the horizon. However, a small 
positive value of AiV ex would imply insufficient amount of inflation after the scale ho crossed the 
horizon. In our numerical work we therefore exclude scenarios with —5 < AiV cx < 50, i.e. we 
(rather conservatively) demand that at least 50 e-folds of inflation can occur after k crossed 
out of the horizon. 

A second consistency condition we impose is that \L\ should not become too large. Specif- 
ically, we require 1-^(^)1 < 20 for all scales between k and the PBH scale. For (much) larger 
values of \L\ our potential (13. 5p may no longer be appropriate, i.e. higher powers of L may 
need to be included. 

Note that eq. ( 13.211) allows to compute the effective spectral index ris{k) exactly: 

m , r ,m u L(k) + l + l\L(k)^ + 3L(k) + l} 

M*) - 1 = M*»0 - 1] Lo + 1 + i . M + 3U + l) ■ < 3 - 24 > 

This in turn allows an exact (numerical) calculation of the spectral index n{k): 

n(k) -l^f" [n, s (k') -l]d]nk'. (3.25) 

In our numerical scans of parameter space we noticed that frequently the exact value for 
n(k) at PBH scales differs significantly from the values predicted by the expansion of eq. (l2.13j) : 
similar statements apply to ns{k). This is not very surprising, given that | ln(koR)\ = 39.1 for 
our value of the pivot scale and Mpbh — 10 15 g. In fact, we noticed that even if as(ko) and 
j3s{ko) are both positive, ns{k) may not grow monotonically with increasing k. In some cases 
ns(k) computed according to eq.( l3.24p even becomes quite large at values of k some 5 or 10 
e-folds below the PBH scale. This is problematic, since our calculation is based on the slow-roll 
approximation, which no longer works if rig — 1 becomes too large. We therefore demanded 
n s(k) < 2 for all scales up to the PBH scale; the first eq. (l3.7p shows that this corresponds to 
i] < 0.5. This last requirement turns out to be the most constraining one when looking for 
combinations of parameters that give large n(fcpBH)- 
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Figure 2: Scatter plot of (3s(ko) vs. as(ko). Here the model parameters c = 2cMj>/Vb, g/c and 
Lq are scanned randomly, with flat probability distribution functions. 

4 Numerical Results 

We are now ready to present some numerical results. We begin in figure El which shows 
a scatter plot of the spectral parameters as(k ) and (3s(k ), which has been obtained by 
randomly choosing model parameters c [defined in eq. fl3.19p ]. g/c and L in the ranged 16 ! 
|c| < 1, \g\Mp/V < 1, \L \ < 20. We require that n s (k ) and as{k ) lie within their 2 a 
ranges, and impose the consistency conditions discussed above. The plot shows a very strong 
correlation between (3s an d as'- if the latter is negative or small, the former is also small in mag- 
nitude. Moreover, there are few points at large as, and even there most allowed combinations 
of parameters lead to very small (3 s- The accumulation of points at small as can be understood 
from our earlier result (I3.13p . which showed that as is naturally of order (ns — l) 2 < 0.004 
within 2 a. Moreover, (3s is naturally of order a 3 J 2 . On the other hand, for as values close to 
the upper end of the current 2 a range, we do find some scenarios where (3s is sufficiently large 
to allow the formation of 10 15 g PBHs. 

We also explored the correlation between ns(fc ) and as(A; ) (not shown). Here the only 
notable feature is the lower bound (13 . 1 6[) on as(k ); values of as(k ) U P to (and well beyond) 

16 We actually only find acceptable solutions for |c| < 0.8. 
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Figure 3: The potential parameters c (in units of Vo/(2M|); double-dotted [green] curve), 
g/(2c) (dotted, blue), the effective spectral index ns — 1 at the PBH scale (dot-dashed, red) and 
the spectral index n— 1 at the PBH scale (solid, black) are shown as functions of Lq = ln(0o/0*), 
for ns(ko) = 0.964 and as(k ) = 0.012. If both solutions in eq.f l3.14p for g/c are acceptable, we 
have taken the one giving larger n at the PBH scale. 

its observational upper bound can be realized in this model for any value of ns(ko) within the 
presently allowed range. Similarly, we do not find any correlation between ns(k ) and fisiko). 
This lack of correlation can be explained through the denominator in eq. fl3.13p . which also 
appears (to the third power) in the expression for (3s once eq. (13.121) has been used to trade c for 
ns — 1: this denominator can be made small through a cancellation, allowing sizable as even 
if ns is very close to 1. Since the same denominator appears (albeit with different power) in 
the expressions for as and fls, it does not destroy the correlation between these two quantities 
discussed in the previous paragraph. 

Figure [2] indicates that large values of n at the PBH scale can be achieved only if as at the 
CMB scale is positive and not too small. In figure [3] we therefore explore the dependence of 
the potential parameters, and of the spectral parameters at the PBH scale on L , for ris(k ) = 
0.964, as(ko) = 0.012. Note that varying L also changes the parameters c and g (or g/c), 
see eqs. fl3.12l) and f|3.14j) . The latter in general has two solutions; however, for most values 
of L , only one of them leads to sufficient inflation while keeping \L\ < 20; if both solutions 
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are allowed, we take the one giving a larger spectral index at the PBH scale, taken to be 
1.5 x 10 15 Mpc -1 corresponding to Mpbh = 10 15 g. Note that ns and n at the PBH scale are 
calculated exactly, using eqs. fl3.7p and f)3.25p . We find that the expansion f )2.13p is frequently 
very unreliable, e.g. giving the wrong sign for n — 1 at the PBH scale for Lq < —1. 

Figure |3] shows that c is usually well below 1, as expected from the fact that n s — 1 oc c, 
see the first eq.( 13.10l) . Moreover, in most of the parameter space eq.f l3.14p implies \g\ < \c\; 
recall that this is also expected, since g is a two-loop term. We find \g\ > \c\ only if |c| is small. 
In particular, the poles in g/(2c) shown in figure [3] occur only where c vanishes; note that the 
spectral parameters remain smooth across these "poles" . 

There are a couple of real discontinuities in figure [3j where the curves switch between the 
two solutions of eq. ( 13. 141) . The first occurs at Lq ~ —0.756. For smaller values of Lq, the 
solution giving the smaller \g/c\ violates our slow-roll condition \ng — 1| < 1 at scales close to 
the PBH scale. For larger Lq this condition is satisfied. Just above the discontinuity, where ns 
is close to 2 at the PBH scale, we find the largest spectral index at the PBH scale, which is 
close to 0.47. Recall from figured] that this will generate sufficiently large density perturbations 
to allow the formation of PBHs with Mpbh = 10 15 g. However, the formation of PBHs with 
this mass is possible only for a narrow range of L , roughly —0.756 < Lq < —0.739. 

At L = —0.31, c goes through zero, giving a pole in g/c as discussed above. Then, at 
Lq = —0.214, the second discontinuity occurs. Here the curves switch between the two solutions 
of eq. fl3.14p simply because the second solution gives a larger spectral index at the PBH scale. 
Right at the discontinuity both solutions give the same spectral index, i.e. the curve depicting 
tt-pbh remains continuous; however, c jumps from about 0.132 to —0.339. The effective spectral 
index ns at the PBH scale also shows a small discontinuity. Recall from our discussion of 
eq.f l2.15p that ns will generally be larger than n at the PBH scale, but the difference between 
the two depends on the model parameters. 

For very small values of |Lo|, c becomes very large; this region of parameter space is therefore 
somewhat pathological. For sizably positive L , n at the PBH scale increases slowly with 
increasing L , while |c| and \g/c\ both decrease. However, the spectral index at the PBH scale 
remains below the critical value for the formation of long-lived PBHs. 

Note that L always maintains its sign during inflation, since L = corresponds to a sta- 
tionary point of the potential, which the (classical) inflaton trajectory cannot cross. For most 
of the parameter space shown in figure El \L\ decreases during inflation. If Lq < decreasing 
\L\ corresponds to V'(<po) < 0, i.e. the inflaton rolls towards a minimum of the potential at 
cf) = 0*. For Lq > we instead have V'(<J)q) > 0, i.e. the inflaton rolls away from a maximum 
of the potential. 

In fact, this latter situation also describes the branch of figure [3] giving the largest spectral 
index at the PBH scale; since here Lq < 0, \L\ increases during inflation on this branch. This 
is illustrated in figures HI which show the (rescaled) inflaton potential as well as the effective 
spectral index as function of either the inflaton field (left frame) or of the scale k (right frame). 
Note that all quantities shown here are dimensionless, and are determined uniquely by the 
dimensionless parameters c defined in eq.f l3.19p . g/c and Lq. This leaves two dimensionful 
quantities undetermined, e.g. Vq and 0*; one combination of these quantities can be fixed 
via the normalization of the CMB power spectrum, leaving one parameter undetermined (and 
irrelevant for our discussion). 

The left frame shows that the (inverse) scale k first increases quickly as <fi rolls down from 
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Figure 4: Evolution of the rescaled inflaton potential, 4Mp(V — Vq) / (Vq<Pq) (solid, black), and of 
the effective spectral index rig (dashed, red), as a function of the inflaton field <f)/<fro (left frame) 
or of the ratio of scales k/ko (right frame). In the left frame the dot-dashed (blue) curve shows 
k/ko, whereas in the right frame it depicts <f)/<f)o] this curve in both cases refers to the scale to 
the right. We took c = —0.1711, g/c = 0.09648, Lq = —0.756; these parameters maximize the 
spectral index at the PBH scale for ris(k ) = 0.964, as(k ) = 0.012 (see figure [3]). 

its initial value 0o- This means that initially moves rather slowly, as can also be seen in the 
right frame. Since as > 0, the effective spectral index increases with increasing k. The right 
frame shows that this evolution is quite nonlinear, although for k/ko < 10 10 , risik) is to good 
approximation a parabolic function of ln(A;/A;o). However, for even smaller scales, i.e. larger k, 
the rate of growth of n$ decreases again, such that n$ — 1 reaches a value very close to 1 at the 
scale k = 1.6 x 10 16 ko relevant for the formation of PBHs with Mpbh — 10 15 g. Recall that we 
only allow solutions where ris(k) < 2 for the entire range of k considered; figures |4] therefore 
illustrate our earlier statement that this constraint limits the size of the spectral index at PBH 
scales. 

The left frame of figure H] shows that the inflaton potential as written becomes unbounded 
from below for — > +oo. This can be cured by introducing a quartic (or higher) term in the 
inflaton potential; the coefficient of this term should be chosen sufficiently small not to affect 
the discussion at the values of of interest to us. Note also that this pathology of our inflaton 
potential is not visible in the right frame, since assuming = 0o < 0* at k = ko, the inflaton 
field can never have been larger than 0*: as noted above, it cannot have moved across the 
maximum of the potential. 

Figure [2] indicated a strong dependence of the maximal spectral index at PBH scales on 
a>s{ko)- This is confirmed by figure 0, which shows the maximal possible n(kp-Qn) consistent 
with our constraints as function of a>s{ko), as well as the corresponding values of the parameters 
c, g/c and L . We saw in the discussion of eq. fl3.14l) that as < is only allowed for a narrow 
range of Lq. In this very constrained corner of parameter space, n(kpBii) remains less than 1, 
although the effective spectral index ns(kpBii) can exceed 1 for as(ko) > —1.7 x 10~ 4 ; recall 
that for the given choice ns(k ) = 0.964, solutions only exist if as(ko) > —3.24 x 10~ 4 , see 
eq. (l3~T6l) . 

For 0:5(^0) < 1.5 x 10 the optimal set of parameters lies well inside the region of parameter 
space delineated by our constraints. ns(kpBn) therefore grows very fast with increasing a^&o). 
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Figure 5: The solid (black) curve shows the maximal spectral index at the PBH scale /cpbh = 
1.6 x 10 19 fco that is consistent with the constraints we impose, as function of as(ko)- The other 
curves show the corresponding model parameters: c (double-dotted, green), multiplied with —5 
for ease of presentation; —g/{2c) (dotted, blue); and — L (dot-dashed, red). ns(k ) is fixed at 
0.964, but the bound on nppH is almost independent of this choice. 

For these very small values of ots(ko), the largest spectral index at PBH scales is always found 
for g = —2c, which implies that the second derivative of the inflaton potential also vanishes at 
= 0*, i.e. 0* corresponds to a saddle point, rather than an extremum, of the potential. 

For slightly larger values of as{ko) the choice g = —2c allows less than 50 e-folds of inflation 
after ko exited the horizon. Requiring at least 50 e-folds of inflation therefore leads to a kink in 
the curve for n(k-pBn)- The optimal allowed parameter set now has considerable smaller \g/c\, 
but larger |c|. 

The curve for the maximal n(fcpBH) shows a second kink at ots{ko) = 0.001. To the right 
of this point the most important constraint is our requirement that \ng — 1| < 1 at all scales 
up to kpBH, as discussed in connection with figures |3] and HI Note that for Q S (fc ) > 0.014, the 
optimal parameter choice leads to % reaching its maximum at some intermediate k close to, 
but smaller than, /cpbh- This leads to a further flattening of the increase of n(kpBK)- 

We nevertheless see that for values of as(ko) close to the upper end of the 2 a range specified 
in ( I3.17P the spectral index at the scale relevant for the formation of 10 15 g PBHs can be well 
above the minimum for PBH formation found in section 2. Figure [U then implies that the 
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Figure 6: The minimal value of ots{ko) that allows the formation of primordial black holes of a 
given mass as a result of density perturbations produced during the slow-roll phase of inflation. 



formation of considerably heavier PBHs might be possible in running mass inflation. However, 
larger PBH masses correspond to smaller /cpbh ; see eq. fl2.10p . This in turn allows for less running 
of the spectral index. In order to check whether even heavier PBHs might be formed during 
the slow-roll phase of running mass inflation, one therefore has to re-optimize the parameters 
for different choices of /cpbh- 

The result of such an analysis is shown in figure [61 We see that in the given model, the 
formation of PBHs that are sufficiently massive, and hence long-lived, to be CDM candidates 
could only have been triggered during the slow-roll phase of inflation if as(ko) is more than 
one standard deviation above its central value. If as(ko) is at the upper end of its present 2 a 
range, PBHs with mass up to 5 x 10 16 g could have formed as result of density perturbations 
created at the end of slow-roll inflation. These results are again insensitive to the value of 
n s (k ). 

5 Summary and Conclusions 

In this paper we have investigated the formation of primordial black holes in the radiation- 
dominated era just after inflation. We have focused on density perturbations originating from 
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the slow-roll phase of inflation. In section 2 we reviewed the Press-Schechter type formalism 
for PBH formation. We found that the formation of PBHs with mass larger than 10 15 g, which 
could form (part of) the cold Dark Matter in the Universe, the spectral index at scale &pbh 
should be at least 1.37, even for the lower value of 1/3 for the threshold <5 t h- This value is higher 
that the value 1.25 found in [27] because here we have assumed that the mass of the collapsed 
region to form PBH is only 20% of the entire energy density inside the particle horizon. 

This spectral index is much above the value measured in the CMB. PBH formation therefore 
requires significant positive running of the spectral index when k is increased. In contrast, 
current observations favor a negative first derivative of the spectral index. For the central value 
n s(ko) — 0.964, where k is the CMB pivot scale, the first derivative as(k ) would need to 
exceed 0.020 if it alone were responsible for the required increase of the spectral index; this is 
outside the current 2 a range for this quantity. However, the second derivative (the "running of 
the running") is currently only very weakly constrained. We showed in a model-independent 
analysis that this easily allows values of n(fcpBH) large enough for PBH formation, even if the 
first derivative of the spectral index is negative at CMB scales. 

In section 3 we applied this formalism to the "running-mass" model of inflation, which had 
previously been shown to allow sufficiently large positive running of the spectral index to yield 
PBH formation. In contrast to previous analyses, we included a term quadratic in the logarithm 
of the field, i.e. we included the "running of the running" of the inflaton mass along with the 
"running of the running" of the spectral index. We showed that this model can accommodate 
a sizably positive second derivative of the spectral index at PBH scales. However, this is only 
possible if the first derivative is also positive and sizable. In fact, like most inflationary scenarios 
with a smooth potential [28] the model does not permit large negative running of the spectral 
index at CMB scales. Current data prefere negative running of the spectral index at CMB 
scales, but are consistent with sufficiently large positive running at the 2 a level. Moreover, we 
saw that a quadratic (in In k) extrapolation of the spectral index to PBH scales is not reliable, 
and therefore computed the spectral index exactly. Imposing several consistency conditions, we 
found that density perturbations that are sufficiently large to trigger PBH formation only occur 
for a very narrow region of parameter space. Among other things, the signs of the parameters 
of the inflaton potential must be chosen such that the potential has a local maximum, and the 
initial value of the field must be slightly (by typically less than one e-fold) below this maximum. 

We emphasize that obtaining sufficiently large density perturbations at small scales is only 
a necessary condition for successful PBH formation. One major challenge of this model is that 
for parameters allowing PBH formation the spectral index keep increasing at yet smaller scales. 
Parameters that lead to the formation of many PBHs with mass around 10 15 g, which could 
form the Dark Matter in the universe, would predict the over-production of unstable PBHs, in 
conflict with data e.g. from the non-observation of black hole evaporation and (for yet smaller 
masses) Big Bang Nucleosynthesis. This problem seems quite generic for this mechanism. One 
way to solve it might be to abruptly cut off inflation just after the scales relevant for the 
formation of the desired PBHs leave the horizon, which could e.g. be achieved by triggering 
the waterfall field in hybrid inflation. However, this is somewhat in conflict with the use of 
the slow-roll formalism of structure formation, which is usually assumed to require a few more 
e-folds of inflation after the scales of interest left the horizon. Moreover, a sharpe end of inflation 
means that the visible universe only inflated by about 45 e-folds; this solves the horizlatness 
problems only if the reheat temperature was rather low [17J. 
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We can therefore not state with confidence that formation of PBHs as Dark Matter candi- 
dates is possible in running-mass inflation; all we can say is that certain necessary conditions 
can be satisfied. Even that is possible only for a very limited range of parameters. This also 
means that constraints from the over-production of PBHs only rule out a small fraction of the 
otherwise allowed parameter space of this model. 
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